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We consider a class of 5-D brane-world solutions with a power-law warp factor a{y) oc y'' , and 
bulk dilaton with profile (j> oilxiy, where y is the proper distance in the extra dimension. This class 
includes the Heterotic M-theory brane-world of Refs. [1,2] and the Randall-Sundrum (RS) model as 
CO . a limiting case. In general, there are two moduli fields y±, corresponding to the "positions" of two 

' branes (which live at the fixed points of an orbifold compactification) . Classically, the moduli are 

' massless, due to a scaling symmetry of the action. However, in the absence of supersymmetry, they 

develop an effective potential at one loop. Local terms proportional to K±, where K± — q/y± is 
(~| ' the local curvature scale at the location of the corresponding brane, are needed in order to remove 

, the divergences in the effective potential. Such terms break the scaling symmetry and hence they 

may act as stabilizers for the moduli. When the branes are very close to each other, the effective 
CNJ ' potential induced by massless bulk fields behaves like V ~ where d is the separation between 

branes. When the branes are widely separated, the potentials for each one of the moduli generically 
develop a "Coleman- Weinberg" -type behaviour of the form a'^{y±)K'±hi{K±/ ^±), where /x± are 
renormalization scales. In the RS case, the bulk geometry is AdS and K± are equal to a constant, 
independent of the position of the branes, so these terms do not contribute to the mass of the moduli. 
However, for generic warp factor, they provide a simple stabilization mechanism. For q > 10, the 
observed hierarchy can be naturally generated by this potential, giving the lightest modulus a mass 
of order m_ < TeV . 

YITP-01-85 



> 



O 



X 



I. INTRODUCTION 



D , Brane-world scenarios, where two or more parallel branes of dimension 3-1-1 are embedded in a "bulk" of larger 
i-C ' dimension have recently been used in order to construct models of considerable phenomenological interest, both in 
particle physics and in cosmology [3,4]. Usually, these models admit an effective 4D description at low energies, where 
the distances between branes are represented by scalar fields (pi{x'^), (^ = 0, ...,3), called moduh. 

An important example of a brane-world is the Randall-Sundrum I scenario (RS) [4], where the bulk is a shce of a 
^ . 5-dimensional anti-de Sitter space (AdS) bounded by two branes of opposite tension. Matter fields may be restricted 
to live on the branes, but in many extensions of the original model some of them are allowed to live in the bulk too 
[5]. In RS there is a single modulus, called the "radion", related to the thickness of the AdS slice. It is given by 
xj}{x^) ~ inpe'xp[—d{x^) / P\, where d is the physical interbrane distance, ^ is the AdS curvature radius and is the 
Planck mass. Although in the 5D theory all fields are assumed to have physical masses comparable to rup, the 4D 
effective mass of fields living in the negative tension brane, given by m ~ rup < >, turns out to be much smaller. 
This is due to a "redshift" effect induced by the AdS geometry between the branes. The exponential dependence of 
the m on < d > can easily explain a large hierarchy between the Planck scale and the electroweak scale without the 
need of invoking very large numbers, which is of course one of the main assets of the RS construction. 

In general, moduli fields are massless at the classical level. The reason is that they are introduced by substituting 
certain integration constants in the classical solutions (e.g. the interbrane distances) by slowly varying fields. Since 
the solutions exist for a continuous range of the integration constants, it follows that the corresponding potential for 
the moduli must be completely flat. Masslessness, however, is phenomenologically undesirable because it gives rise 
to long range scalar interactions which are severely constrained by observations, and therefore it becomes necessary 
to introduce a stabilization mechanism. Needless to say, this can be done in an ad hoc manner, by introducing new 
interactions determining the equilibrium distances between the branes. Alternatively, it is possible that a suitable 
mass may be generated by quantum effects. 

In the absence of supersymmetry, moduli fields tend to develop an effective potential. Of course, this happens 
already in the simplest Kaluza-Klein (KK) compactification, and even if there are no branes [6]. A five dimensional 
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field X can be split in an infinite tower of massive KK fields, labeled by a discrete index n. The masses of the KK 
fields m„((/3) depend on the size <p of the extra dimension, and because of that an effective potential V(ip) is generated 
at one loop. From the 5D perspective, this corresponds to the (nonlocal) Casimir energy density due to the presence 
of compact directions. 

In most cases, the self-gravity of brane and bulk matter content induce a warp in the extra dimension. In this 
situation the subtraction of infinities in the calculation of V{ipi) is more elaborate than it is in flat space. In particular, 
local terms proportional to worldvolume operators on the brane are generated, which may give rise to interesting (and 
sometimes surprising) physical effects. This paper follows up on our previous work [7], where we considered the 
particular case of the Randall-Simdrum (RS) model. A curious result [7] was the absence of logarithmic terms in the 
effective potential. As mentioned above, in the 4D effective description, the fields living in the negative tension brane 
have masses which are proportional to the radion expectation value. This is reminiscent of the Higgs mechanism 
in the standard model, and naively one might have expected the usual 4D logarithmic effective potential for the 
radion. However, we showed that a regularization which would preserve the five dimensional general covariance did 
not produce such logarithm. Consequently [7], it was found that the one loop effective potential is not sufficient to 
stabilize the radion at a phenomenologically acceptable mass in this model. Several papers have later appeared in the 
literature [8-11] where the effective potential is obtained in dimensional regularization, through certain subtractions 
corresponding to a renormalization of the brane tensions in the dimensionally extended spacetime. The results 
obtained in dimensional regularization were eventually found to coincide with ours, giving an indication that both 
procedures are in fact equivalent. The RS model, however, is somewhat special in that the bulk and branes are 
maximally symmetric, and all possible coimterterms amount to renormalizations of the brane tensions. Therefore, 
in order to unveil the generic properties of the effective potential in the presence of self gravitating branes, the 
consideration of more general cases seems to be needed . 

In this paper we shall study a more general class of warped Brane- World models, in which the bulk is no longer 
AdS. The model we consider contains a bulk "dilaton" scalar field with an exponential potential which is coupled 
to the 5-dimensional Einstein gravity. This provides a family of solutions with two classically massless moduli (a 
combination of which is the distance between branes). The solutions have a power- law warp factor a{y) oc in place 
of the exponential warp of the RS model. Here y is the proper distance in the extra dimension and g is a constant 
which depends on the parameters in the Lagrangian. For g = 1/6 this reduces to the Heterotic M-theory brane- world 
of Ref. [1], which may perhaps be relevant for the recently proposed Ekpyrotic universe scenario [2], whereas the RS 
model is recovered in the limit oo. 

The plan of the paper is the following. Section II is devoted to the classical dynamics of the model. After 
introducing the background solution, in Section II. b we derive the action for the 4D moduli fields. These are massless 
at the classical level, due to a scaling symmetry of the 5D action which we discuss in Section II. b. The hierarchy 
between the Planck scale and electroweak scale is discussed in Section II. c, since in general there are some differences 
between the present case and the RS case. In Section HI we give a brief account of quantum effects by considering 
the effective potential induced by a conformally invariant field. This case is rather trivial, and it properly illustrates 
the Casimir interaction between the branes. However, it misses the possibility of interesting local terms induced at 
one loop. Section IV deals with the formalism for the calculation of V{ifi) in more general cases, and contains the 
core of the discussion on renormalization. We start, in Section IV. A, with a digression on the specification of the 
measure of functional integration in a curved background. General covariance alone is not sufficient to determine the 
measure in the present situation, due to the presence of a nontrivial dilaton. Because of that, we shall consider a 
one parameter family of possibilities, related with each other through dilaton dependent conformal transformations. 
Sections IV. B, C and D are devoted to dimensional regularization and zeta function regularization, which are shown 
to give equivalent results. The explicit calculation of the effective potential in some simple cases is given in Section 
V. The Heterotic M-theory model is briefly discussed in Section VI. Possible consequences for the stabilization of the 
moduli are discussed in Section VII, and our conclusions are summarized in Section VIII. Some technical issues are 
left to the Appendices. 



The classical backgroimd we shall consider in this paper is a 5-dimensional spacetime with a nontrivial background 
scalar field (/), which we shall refer to as the "dilaton" . The fifth dimension is compactified on a Z2 orbifold with two 
branes at the fixed points of the Z2 symmetry. The action for the background fields is given by 
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(7+ J d'^x^e"^''^ - a- jd^x^- 



e^^/^ (1) 



where TZ is the curvature scalar, K5 = IGttGs, where G5 is the 5-dimensional gravitational coupling constant. We 
have denoted the induced metrics on the positive and negative tension branes by and g~^,, respectively. To find a 
solution of the equations of motion, we make an ansatz where the 4-dimensional metric is flat, 

rfs^ = dy"^ + a^{y)'n^^dx^'dx", (2) 

with a a;''-independent scalar field (j) = <p{y). The positive and negative branes are placed at y = y+ and y_, 
respectively. Under these assumptions, the equations of motion for (a, (f)) in the bulk become 

4> + A^^=U'{(t>), (3) 

where U{(j)) = Ae"^"^, a dot represents differentiation with respect to y and a prime represents differentiation with 
respect to (p. 

As shown in [13], there is a solution of Eqs. (3) for any value of c given by 



where 



^Gq ln(y/?/o), 

a{y) = {y/yor, (4) 



2 / 3g(l-4g) 

(Constant rescalings of the warp factor arc of course allowed, but unless otherwise stated, we shall take the convention 
that a{y) = 1 at y = yo.) Assuming y- < the boundary conditions which follow from Z2-symmetry imposed on 
both branes are given by 

0± = T^^±e(^/^)*±, (6) 
^=±iK5a±e('=/2)*±, (7) 



6^ 
a 



and they are satisfied if a± are tuned to 



, 1 / 48gA 
K5 V 1 - 4g 



In the absence of the branes, the spacetime (4) contains a singularity at j/ = 0. Since we are considering the range 
between and y+, this singularity is of course inocuous. Our spacetime consists of two copies of the slice comprised 
between y_ and y+, which are glued together at the branes. Hence, the 5-th dimension is topologically an S^/Z2 
orbifold. 

For (7 = 1/6 this solution is precisely the M-theory heterotic brane of Ref. [1] . On the other hand, the RS case, where 
the bulk is AdS and there is no scalar field, can be obtained by taking the limit g — > 00 and j/o — > 00 simultaneously, 
while its ratio is kept fixed, 



hm ^ = J^. (9) 

g^oo g V A 



Defining y = yo + y* , we find that in the limit the warp factor becomes an exponential 

lim a = e^'l\ 

q—>oo 

which corresponds to AdS space with curvature radius equal to £. 
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A. Moduli fields 



For fixed value of the coupling c, the solution given above contains only two physically meaningful free parameters, 
which arc the locations of the brancs y_, and This leads to the existence of the corresponding moduli, which 
are massless scalar fields from the 4-dimensional point of view. In addition to these moduli, the massless sector also 
contains the graviton zero mode. To account for it, we generalize our metric ansatz (2) by promoting ry^y to an 
arbitrary four-dimensional metric: 

ds^ = dy^ + a'^{y)g^,4x)dx^dx''. (10) 

It is easy to show that TZ = a~'^iZ + TZ^^^ , where TZ^°^ is the background Ricci scalar in five dimensions and TZ is the 
four-dimensional one. For constant values of the metric and moduli, we have a solution of the equations of motion 
whose action vanishes. Hence, only the terms which depend on derivatives of the metric or derivatives of the moduli 
will survive after the five-dimensional integration. This fact can be used in order to simplify the derivation of the 
effective action for the moduli, since all terms without any derivatives can be dropped. 

In the bulk, all terms will cancel except for the one which is proportional to the four- dimensional Ricci scalar, "R. 
Let us consider the contribution from the branes. The metric induced on the branes is given by 

a'^u = a±[ff;*L/ + a±^d^,y±d^y±\. 

Here, and in what follows, the subindices ± mean that the quantity is evaluated at the perturbed brane location. The 
brane tension terms in the action contain the determinant 



V-9± = a±V^ 
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which induces kinetic terms for the moduli. Here, the tilde on the kinetic term indicates that the derivatives are 
contracted with the metric g. In fact, the five dimensional Ricci tensor TZ contains second derivatives of the metric 
and therefore it is singular on the brane, giving a finite contribution to the action. To handle this contribution it 
is simplest to introduce fiducial boundaries in the neighborhood of the branes, were we add (back to back) pairs of 
Gibbons-Hawking boundary terms. These have the form 



- / d^x 

K5 J 



9^IC±, (11) 



where /C± is the trace of the extrinsic curvature of the fiducial boundary near each one of the branes. The action is 
separated into two parts. The first one is a "bulk" part, consisting of an integral over two copies of an open set which 
excludes the branes, (y_ -I- e) < y < {y+ — e), plus terms of the form (11) at the boundaries y = y- + e and y = y+ — e. 
Then there is a "brane" part, which includes an integral of the action over the infinitessimal open sets of thickness 
e around the branes, suplemented with terms of the form (11) at the boundaries of these open sets (these "brane" 
boundary terms sit back to back with the ones used in order to bound the "bulk" , and have opposite sign relative to 
them, since the normal to the fiducial boundary has opposite sign on each side of the boundary. Thus the total effect 
of the boundary terms is to add zero to the action). As is well known, through integration by parts the boundary 
terms absorb the second derivatives in the Einstein term TZ. Hence, the singular contribution from the gravity kinetic 
term on the brane disappears, and in the limit e ^ 0, the only contribution to the action from the interval of width 
2e around the branes (together with the added boundary terms), is from the brane tension term itself, but not from 
the gravity kinetic term. On the other hand, the bulk contribution has to be supplemented with the boundary terms 
(11). 

The boundary terms can be evaluated as follows. Consider, for instance the hypersurface which is located at 
y = y+{x'^). In terms of the new coordinate y = y — y+{x'^), the brane is at ^ = 0, and the metric can be written as 

ds'^ = N'^df + g+^iN^'dy + dx^'){N''dy + dx"), (12) 

where we have introduced the lapse function N'^ = 1 — g'^^N^N", and the shift vector = g^'^'^y^v Then, the trace 
of the extrinsic curvature is given by 



/ 



= / dx' 



/■ 



2N 
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(13) 



where the vertical hne means a covariant differentiation with respect to the induced metric 5+,. We neglected the 
terms which are higher order in in derivatives of the modulus y+. As mentioned before, the subindices ± mean that 
the quantities arc evaluated at the actual position of the brane, and the expression is in fact nonperturbative in the 
positions y± themselves (although not in the derivatives). 

Substituting the previous expressions into the action (1), with the addition of the extrinsic curvature terms, and 
using the background equations of motion we find 



n 



This can be rewritten as 
Here we have introduced 



-1 



.2/0 



6q 



9 + 1/2 



9+1/2 



(14) 



(15) 



and the four dimensional Newton's constant G given by 



G=[q 



yo 



(16) 



The modulus corresponding to the positive tension brane has a kinetic term with the "wrong" sign. However, this 
does not necessarily signal an instability, because it is written in a Brans-Dicke frame. One may go to the Einstein 
frame by a conformal transformation. It is convenient to introduce the new moduli y and ij) through [2] 



ip+ = ip coshtp, (p- = ipsinhip, 



and to define the new metric 



It is then straightforward to show that ^/g TZ = \f—g + 6 <^ ^ (^v)^]- Substituting into the background action 

(15), we have 
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1 + 2gr 



l + 2g 



(17) 



Therefore, both moduli have positive kinetic terms in the Einstein frame. At the classical level, the moduli are 
massless, but as we shall see in the following Sections, a potential term of the form 



(18) 



is generated at one loop, which should be added to (17). 

In the RS limit q ^ 00 [see Eq. (9)] the kinetic term for one of the moduli disappears. This is to be expected, 
because the bulk is the maximally symmetric AdS space. In this case only the relative position of the branes y+ — y- 
is physically meaningful and the other modulus can be gauged away (see also [12] for a recent discussion of this case). 

Matter does not couple universally to the 4D Einstein metric g^^. Rather, it will couple to combinations of the 
metric and the moduli fields. For completeness, in Appendix A we give an alternative form for the moduli action in 
terms of the induced metrics on either brane, to which brane matter fields are universally coupled. 
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B. Scaling symmetry 



Since we are interested in the effective potential for the moduH, it is perhaps pertinent to start by asking why these 
fields are massless at the classical level. The reason is that under the global transformation 

9ab ^ T^Qab, (19) 



0_^_(2/c)lnr, (20) 

the action (1) scales by a constant factor 

Here gat is the metric appearing in the action (1). Acting on a sohition with one branc, the transformation simply 
moves the brane to a different location. Hence, all brane locations are allowed, from which the masslessness of the 
moduli follows. However, we should hasten to add that this is just a global scaling symmetry which need not survive 
quantum corrections. 

It is interesting to observe that by means of a conformal transformation, we may construct a new metric g^^^ which 
is invariant under the scaling symmetry 

9l:,!=e^*9ab. (21) 
In terms of this new metric the action takes the form 

-a+ j d'^x^ - a_ jd^x^ e-^^^^/^. (22) 

Now, the symmetry is a more shift in (f). Moreover, with our background solutions for gab and (j), the metric g'^j^ is 
just AdS, as can be easily shown from (21) and (4). 

For certain discrete values of c, the action (22) can be obtained from dimensional reduction of (5 + n) dimensional 
pure gravity with a cosmological term A, where the additional n dimensions are toroidal [24]. In this case, the factor 
g-3c0/2 ^jjg overall scale of the internal n-dimensional volume, and the value of c is given by = 2/3q, with 

n + 3 



It is not surprising, then, that the metric g^^' in (22) corresponds to the "external" components of a 5 + n dimensional 
anti de Sitter space, since the starting point is in fact pure gravity in (5 + n) dimensions with a negative cosmological 
term. The calculation of quantum corrections in this higher dimensional space, and its relation with the calculation 
of quantum corrections in the effective 5D theory which we consider in this paper, will be reported in a separate 
publication [24]. 



C. The hierarchy 

As mentioned in the introduction, one of the motivations for studying brane-world scenarios has been the search 
for a geometric origin of the hierarchy between the effective Planck scale rup and the electroweak scale. In the 5D 

— 1/3 

description, all matter fields are assumed to have masses which are close to the cut-off scale of the theory M = 
And yet, with the help of an exponential warp factor (as in the RS model) it is easy to generate a hierarchy of the 
order of mp/m ^ 10^^. Here m is the effective mass of fields which live on the negative tension brane, as "seen" in 
the effective four dimensional description [4] . In this subsection we shall review this mechanism, including the case of 
a warp factor with arbitrary power q, since there are some minor differences with the RS case. 
The effective four-dimensional Planck mass is given by 

(23) 



1-1^ 

y+ 



2q+l 
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where M is the 5-dimcnsional Planck mass, as can be seen from Eqs. (15) and (16)], where, without loss of generality, 
we have taken = yo- Here, and for the rest of this Section, we shall follow standard practice and refer all physical 
quantities to the measurements of clocks and rods located on the positive tension branc. 

Let us now consider the mass scales of fields which live on the branes. We expect these fields to couple not only 
to the metric, but also to the background dilaton cj). There are many possible forms for this coupling, but it seems 
reasonable to restrict attention to those which respect the scaling symmetry (19-20). For a free scalar field ^' which 
lives on the negative tension brane, and whose mass parameter is comparable to the cut-off scale, the action takes the 
form 

Here we have introduced a fudge factor a to allow for an intrinsic mass which is slightly lower than the cut-off scale. 
The function F{(j)) can be reabsorbed in a redefinition of \E', and thus the relevant warp which determines the hierarchy 
between mass scales on the positive and in the negative tension branes is the one corresponding to the metric g^*' . 
Our field will be perceived from the point of view of the 4D effective theory as having a mass squared of order 



~ olM^ ( . (25) 



Notice that there arc two different factors which determine the hierarchy between m and m^. The first one is the warp 
factor (y_/y+)5~^ appearing in Eq. (25), which "redshifts" the mass scales of particles on the negative tension brane 
(except for g < 1, in which case the particles on the negative tension brane appear to be heavier than those on the 
positive tension brane). This generates the hierarchy in the RS model. The second one is the possibly large volume 
of the internal space, which enhances the effective Planck scale with respect to the cut-off scale [see Eq. (23)]. This 
generates the hierarchy in the ADD model [3] with large extra dimensions. Considering both effects, the hierarchy h 
is given by 



2q-2 

ml ~ 2 My+ \y+J 

It is known that without a warp factor, it is not possible to generate the desired hierarchy from a single extra 
dimension, since its size would have to be astronomical. An interesting question is what is the minimum value of the 
exponent q which would be sufficient in order to generate a ratio m/nip 10~^^. The best we can do is to take the 
curvature scale {y+/q) slightly below the millimeter scale, 

{y+/q) < mp{TeV)-^ ~ mm. (26) 

in order to pass the short distance tests on deviations from Newton's law. On the other hand, we also need 

{y./q) > M-\ (27) 

since for smaller values of y_ the curvature becomes comparable to the cutoff scale M and the theory cannot be 
trusted ^ Substituting in (23) we have > mp{TeVf and {y-/y+) > {mp/TeV)-'^/^ , which leads to 



W-"~'-^>a\ — ] . (28) 



4(2<;-l) 
3 



V 



Hence, a warp factor with exponent g > 5/4 may account for the observed hierarchy with a single extra dimension, 
but it appears that this cannot be done for lower values of q. ^ In particular, the Heterotic M-theory model, with 
q= 1/6, does not seem to allow for such possibility. 



^Onc should also bear in mind that the scale M might itself be a "derived" quantity, as it happens for instance with the 
Planck mass in theories with additional large extra dimensions. In this case, the true cut-off scale may well be below M. Hence, 
the lower bound (27) on y- should just be considered a necessary condition for the low energy description not to break down. 

^Except, of course, by giving up the assumption that the Lagrangian of matter on the branes should scale in the same way as 
the rest of the classical action [see the discussion around Eq. (24)]. If we allow any coupling of to the mass term for $, then 
any hierarchy can be easily generated for any value of q. 



7 



III. EFFECTIVE POTENTIAL INDUCED BY BULK CONFORMAL FIELDS 



Before embarking on a detailed discussion of renormalization, we shall consider in this Section the case of a conformal 
bulk field X- This case is rather easy to handle, and it is useful in illustrating the non-local contribution of the 
vacuum energy to the masses of the moduli. Unfortunately, it misses possibly interesting effects which may arise due 
to renormalization of local operators on the brane, whose discussion we defer to the next sections. 

Following the discussion given in § III of [7], we define the conformal coordinates by 



dy 



X 1-9 

yo y ' 



(29) 



and we rewrite the metric as 

ds'^ = a^{z){dz'^ + rt^^dx'^dx''), a{z) = [z/zof, (30) 

where 

P = -. > ^0 = -pi r- (31) 

1-9 |1-9| 

Here we should mention that the direction of increasing z does not coincide with the direction of increasing y when 
q>l. 

The Casimir energy density p of the conformally coupled scalar field in this spacetime is related to its counterpart 
in a flat spacetime with a compactified extra dimension (a = 1 in the above metric) by a conformal transformation. 
The relation is p = a~^po, where po is the flat space value: 

""'A^r^A^^' ^-i^<iM)=»2.46.10- (32, 

Here z+ and z_ are the position of the positive and the negative tension branes in the conformal coordinate z. The 
double signs in (32) refer to bosons or fermions respectively. Then, we find that the contribution of a conformally 
coupled scalar field to the effective potential per unit co-moving volume is given by 

V{z+,z.) = ^gV = 2 j a\z)pdz = TT-^—Ti. (33) 

Here, V is the effective potential per unit "physical" volume (i.e. the volume as measured with the Einstein metric 
g), to be inserted in (17). In terms of the moduli fields 93 and ip we have 

V{ip,i)) =T B AV"^^^ [(coshV)^^-^ - (sinhV)^^-^]"\ (34) 

where we have defined 

^^4(1-^ 1 



V(l-4g)2' ' l + 2g' 

and we have used \/g = (^^ for the background solution. It should be mentioned that local terms induced by quantum 
corrections may be added to (34), both due to fields which live on the branes as well as from non-conformal bulk 
fields. A discussion of these terms is deferred to the next section. 

In the RS case 7^0 the ip dependence in (34) disappears, and we recover the results of [7] . We have 

For ip <^1, the field ip corresponds to the hierarchy between scales on both branes: 

, a_ TeV 

tanh^ ip = — ~ <c 1, 

a+ mpi 

where mpi ^ lO^^GeV is the Planck mass. In the second equation we assume that the warp is responsible for the 
hierarchy between Planck and electroweak scales. In this case we have 

Vrs{^Ij) « 4 V^[l + 4^ + ...] ~ A {TeV)\ 
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Here we have assumed that the AdS radius is not far below the Planck scale, in which case the contribution to the 
effective potential is of electroweak order. This huge contribution must somehow be cancelled by some other term in 
order to have an acceptable four dimensional cosmological constant. Also, the slope of the potential has to be cancelled 
at the place where the radion sits. These two conditions can be imposed if we allow for a finite renormalization of 
the brane tensions, which contribute proportionally to oc cosh'' w 1 and al oc sinh'* ip ip"^ for the positive and 
negative tension branes respectively. With these additions we have 

Vrs(V') = At-^ [ci + C2V^ + + ...] , 

where Ci and C2 are undetermined constants which can only be fixed by experiment. The "renormalization" condition 
dti,V{ip) = at the observed value of the radion ip = tpobs ^ TeV/mpi forces C2 ~ V'obs- Hence the mass of the radion 
induced by this quantum correction is given by [7] 

ml^ A -|V - ^ i>obs{TeVf. (35) 
t nipi 

Again, in the last equality, we have taken I to be similar to the Planck scale. The radion mass is much lower than the 
electroweak scale and hence it is not enough for a phenomenologically viable stabilization of the radion. The reason 
for such as small mass is that V(^) is polynomial for small V, with the first power being tp'^. The situation would 
be different if the potential took the form xj}'^ log(V') as in the Coleman- Weinberg case. Then, the induced mass 
would be of the same order as the induced potential. Remarkably, the absence of logarithmic terms seems to extend 
to the contribution from massive fields too [8], and so in the RS case it appears that the radion cannot acquire a 
phenomenologically interesting mass solely from one loop quantum effects. 

For the more general case of a power law warp factor, the possibility of an efficient stabilization by quantum effects 
will depend on the nature of the local operators which are induced by quantum corrections. We now turn to a 
discussion of this subject. 

IV. EFFECTIVE POTENTIAL IN THE GENERAL CASE 

In this Section we set up the framework for computing the contribution to the 1-loop effective potential from a scalar 
field X propagating in the bulk with a generic mass term, which may include couplings to the curvature of spacetime 
as well as couplings to the background dilaton cf). The effective potential for the moduli y± will be defined as usual in 
terms of a Gaussian path integral around the background solution. Before presenting the actual calculation, however, 
a digression on the choice of the measure of integration will be useful. 

A. Specification of the functional mesisure 

A quantum field theory is defined not just by the classical Lagrangian, but it is also necessary to specify the 
measure of functional integration. The latter is usually prescribed by demanding certain symmetries or invariances. 
For instance, for scalar fields in curved space, invariance under diffeomorphisms is an obvious requirement. If gravity 
were the only background field, then this requirement would suffice to uniquely define Gaussian integration around 
that background. On the other hand, if there are fields other than gravity with a nontrivial profile (such as our dilaton 
(/)), then there is a wide class of possible choices, related to each other by dilaton dependent conformal transformations. 
All choices within this class are equally good from the point of view of diffeomorphism invariance. Strictly speaking, 
however, they are inequivalent due to the well known conformal anomaly. 

To be definite, let us concentrate in the simple case of a bulk scalar field x with canonical kinetic term. The 
(Euclidean) action for this field is given by 

S[x]^\jd^x^gxPx, (36) 
where we have introduced the covariant operator 

P = -{Ug+E). 

Here Dg is the d'Alembertian operator associated with the metric g^b, and E = -E[9a6i'/'] is a generic "mass" term. 
Typically, this takes the form E = —m^ — £,Tlg, where m is a constant mass, TZg is the curvature scalar and ^ is an 
arbitrary coupling. Throughout this Section we shall leave E unspecified. 



9 



A volume measure Vx in field space can be found from a metric Gxy on J^, through the relation: 

Vx = ^/Gl[dx''. (37) 

X 

Here, the spacetimc coordinates x and y are considered as continuous labels for the coordinates = xi^) of the 
infinite dimensional space JF, and G is the determinant of G^y To specify G^y, we note that a natural definition of a 
scalar product in the space of field variations 6x can be given in terms of the spacetime measure dfjb{x), through the 
relation 



dn{x)dn{y) G^y Sxi Sxl = j diJ,{x)Sxi{x)dx2{x). 



We denote field variations by 5x just to emphasize that we are refering to elements of the tangent space. More precisely, 
^X = Sx'^Sx, where = d/dx^ is the coordinate basis of the tangent space at the point p which corresponds to the 
background solution. In a Riemannian spacetime, the invariant measure is given by 



dii{x) = ^g{x)d°x, (38) 

where g is the determinant of Qah, and D is the dimension. The implicit definition of G^y given above is just the 
identity 5^{x,y) with respect to dfi integration, 



(5(")(a;-y) 
V5(a 



Gxy = 5^{x,y) = "-^^=^. (39) 



It is convenient to express the field variations in an orthonormal basis Xn, with < Xn,Xm >= Snm, so that 5x{x) = 
X^n '^"Xn(^)- this basis, the components of the field variation are c", and the metric is just the usual delta function 
(the continuous or the discrete delta function depending on whether the normalization of Xn is continuous or discrete): 

Gnm — ^nm- 

Substituting in (37), we have 

Vx = Jldc^. (40) 

n 

It should be clear from the previous discussion that the definition of T^x is associated with a natural definition of 
d^{x). However, in the problem under consideration in this paper, the choice of d/U is not unique. In our case, there 
is a nontrivial dilaton field (j), and we can consider a whole class of spacetime measures of the form 

d^ie{x) = y/g^ d^x = ng{4>)^ d^x, 

which correspond to conformally related metrics 

for an arbitrary fimction ilei<f>) (we shall use some specific choices of this function in the next subsection). In the 
presence of a dilaton, the coupling to gravity is not universal and it is not clear which one of these metrics should be 
considered more physical. To proceed, it is convenient to define the operator Pe associated with the metric g^f^ by 

= n,-2p, (41) 

where P was introduced after Eq. (36). This operator can be written in covariant form as 

Pe = -{ae + Ee), 

where 

Ee = (^^^) °elnl^e - {^Y^) < da\^^e dblnQe + nfE, 
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and Og is the covariant d'Alembertian corresponding to g^^,. Introducing xe = ^^^^X> the action for the scalar 
field can be expressed as 

^[X] = IJ d^x^e xe Pe Xe- (42) 

In terms of the field xe has a perfectly canonical and covariant kinetic term. 

Thus, the same arguments which lead to (39) can now be used in order to find the natural line element in field 
space associated with the spacetime measure dfjb0{x): 

dSg = j j diieix)dfxeiy)Gly dxe dx0 + ■■■ = j d'^x^ge^x) {dxe{x)f + ... 

Here, the ellipsis denote the omitted terms which correspond to variations of other fields in the theory (in particular, 
these include the variations of the gravitational field and the dilaton) . Let us compare this line element with the one 
considered above 



dS^^ j j d^{x)d^i{y)G,y dx^ dx^ + ... = j d^xyMxj {dx{x)Y + ... 

For field variations where x changes but the rest of the fields (metric, dilaton, etc.) are constant, we have 

dSl = j d^'x nl^/^ {dx{x) f + dS" (5, ... = cansU dx^ + 0), (43) 

and therefore dS^ 7^ dS^ in general. Of course, the corresponding measures of integration will also be different. In 
the basis {x0n} which is orthonormal with respect d/Ug, the field variation can be expanded as (5xe(a;) = c^XSm 
and the new measure takes the form 

^T^x)e=Yi'^c^o- (44) 

n 

Using xem = ^e^^'^Xm, it is straigtforward to show that c™ = M™Cg, where M™ =< Xm,^0^Xn >^i= i^e^)7^- 
Hence the two measures (40) and (44) are related by 

2?X = Je {Vx)e, (45) 

where the Jacobian is formally given by 

J0 = deting^) = exp [-TrlnOe] . (46) 
In the last step we have used the formal definition of the L2 trace: ^ : 

= J2[d''x XmiOxm) = E / ^""^ 9e'^^ Xem{Oxem). 

m •' m 

The trace is well defined if the diagonal matrix elements of the operator O decay sufficiently fast at large momenta. 
Unfortunately, the diagonal matrix elements of In fig do not decay at all at large m, and so the trace is ill defined 
unless we introduce a regulator. We will address this question below, where we will explicitly define what we mean 
by Je- 

Perhaps we should add, for the sake of clarity, that the difference between the line elements dS"^ and dSg, and 
consequently the difference between the associated measures, is not due to field redefinitions. Both objects are 
different, but since they are defined geometrically, they are both invariant under field redefinitions (in the same 
sense that any line element is invariant under coordinate transformations). Rather, the relation (46) expresses the 



^The definition of the trace is robust, in the sense that it is independent on the metric one uses in order to define the 
orthonormal basis, as long as the corresponding measures are in the same L2 class. This will be the case, for instance, if the 
metrics are related by a conformal factor which is bounded above and below on the manifold. 
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well known conformal anomaly. The measure is not invariant under conformal transformations because these do not 
correspond to a change of coordinates in field space T . They correspond to a change of the spacetime metric and 
consequently to a change of the metric on T . This sort of ambiguity does not arise when we consider scalar fields in 
flat space. Consider an action with a general kinetic term of the form: 



A natural line element in field space can be obtained by "stripping off" the flat metric if^ and replacing the partial 
derivatives with differentials of the fields: 



This procedure cannot be transported into a curved space, because in erasing the factor g"*" from the kinetic term, 
it makes a difference what exactly we have chosen to call the metric of spacetime: gab or (this is, by the way, 
the reason why the factor of fig appears in (43)). Usually, flat space definitions can be generalized to curved space 
through the principle of general covariancc: objects should bo defined geometrically, and they should reduce to their 
flat space definition when the spacetime metric is flat. The question is, however, which object should be considered 
to play the role of spacetime metric, so that we know when to call it flat. Physically, too, one should expect that a 
preferred spacetime metric should play a role in regularizing and renormalizing the theory. Suppose that we attempt to 
regularize with a physical cut-off, so that all degrees of freedom beyond a certain scale are ignored. In our background 
(which is conformally flat) a constant physical cut-off scale corresponds to a different co-moving scales at different 
places. The relation between physical and co-moving scale is of course given by the metric, and therefore it can make 
a difference which one we use. 

(s) 

Since wc have a classical scaling symmetry in the gravity and dilaton sector, one could argue that .g^jf, , which 
is invariant under scaling (see Section II. b), is the preferred physical metric. ^ However, even in this case the 
divergent part of the effective potential will not respect the scaling symmetry, and consequently we need to introduce 
countcrtcrms with the "wrong" scaling behaviour. Hence, in what follows, we shall take the conservative attitude 
that the measure is determined in the context of a more fundamental theory (from which our 5-D effective action is 
derived), and we shall formally consider on equal footing all choices associated with metrics in the conformal class of 

(s) 

gab, including of course As we shall see, the difference between these choices amounts to the addition of local 
terms in the effective potential. 

The contribution of the field x to the renormalized effective potential Ve per unit co-moving volume parallel to the 
branes is given by: 

exp [-AiVe + Vfn] ^ J Px)« e-^^^^ = (detP,)-l/^ (47) 

where A is the co-moving volume under consideration and we have used (42) and the measure (44) to express the 
gaussian integral as a determinant. The term V/'" is a local counterterm which, in dimensional regularization, needs 
to be subtracted from the regularized effective potential (its explicit form will be given in the coming sections). In 
zeta function regularization, the left hand side is already finite, and V/™ is unnecessary (it corresponds to a finite 
renormalization of couplings). Eq. (47) can be written as 

F, = ^ln(detP,)-y/-. (48) 



Eq. (45) suggests the notation 



^0 = ^ ln(det P) - Vt + \^^Je- (49) 



The reader should be aware, however, that the definition of the Jacobian in Eq. (46) is only formal because the trace 
in the r.h.s. of this equation is ill defined. For that reason, it is not clear that Eq. (49) would hold with the definition 



"^Notc, in particular, that the overall scaling factor of the action (1) under (19) and (20) depends on the spacetime dimension, 
and hence the symmetry itself is different when we change the dimension. By contrast, the scaling of (22) remains the same in 
any dimension. 
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(46), after substituting determinants by traces and applying any kind of regularization to the formally divergent 
traces. To avoid misinterpretations, in the discussions that follow we shall take Jo to be defined by Eq. (49), that is 




(50) 



where the expression in the right hand side is to be calculated in some regularization scheme. 

The way the 6 dependence of Ve arises is very different in different regularization schemes. In Eq. (49), the 
determinant of P is independent of 9 (we recall that this operator corresponds to the choice = 1). In dimensional 
regularization, In Jg vanishes, but the divergent term Vg'^'" which is subtracted from ln(det P) depends on the choice of 
physical metric g^^. On the other hand, in zeta ftmction regularization, ln(det P) is finite and does not play a role 
(in any case, any finite renormalization does not introduce a dependence in 0). Rather, in this case, the dependence 
on 9 comes from In Jg, which does not vanish in this regularization scheme. In both cases, the 9 dependence of Ve is 
the same. 

As we shall see, this dependence can be cast in the form of local operators on the branes, and therefore the ambiguity 
in the choice of the integration measure can also be understood as modification of the classical action. It should be 
noted, however, that the local operators which result from a shift in 9 have different form than the terms arising 
from the usual shift in the renormalization constant /z which inevitably crops up in the regularized traces. In the 
cases we shall consider, the latter will take the form -ft'^(y±), where K denotes terms which behave like the extrinsic 
curvature of the branes at the positions y±. On the other hand, the 6'-dependent terms behave as K^{y±)<t){y±)- 
Since K{y) behaves like the inverse of y whereas ^{y) behaves logarithmically with y, these terms will give rise to 
Coleman- Weinberg type potentials for the moduli. 



The direct evaluation of the determinant of Pg appearing in Eq. (47) turns out to be rather impractical, due to 
the complicated form of the implicit equation which defines its eigenvalues. For actual calculations it is convenient to 
work with a conformally related operator Pg whose eigenvalues will be related to the KK masses. 

Following [7] , we introduce a one-parameter family of metrics which interpolate between a fictitious flat spacetime 
and any of the metrics in the conformal class of the Einstein metric: 



where 9 parametrizes a path in the spacx; of conformal factors. For definiteness we shall restrict attention to conformal 
factors ((/>) which have an exponential dependence on the dilaton: 



With this choice, 6* = represents flat space and 6* = 1 corresponds to the Einstein frame metric (30). For 9 = 
the metric coincides with the metric g'^J introduced in Subsection II. B, which is invariant under the scaling 
transformation (as mentioned before, this metric corresponds to a five dimensional AdS space, with curvature radius 
given by ^o)- 

The operator Pq = P0=o is the wave operator for the KK modes which one would use in a four-dimensional 
description. The Lorentzian equation of motion Px = can be written as 



B. Conformal transformations and the KK spectrum 




'•ei^) 9ab, 



(51) 




(52) 



PoXo = 



where 



Po = -aD-i + M^z). 



(53) 



Here Dd-i is the flat space d'Alembertian along the branes, and 



M2 = _92 _ 



is the Schrodinger operator whose eigenvalues are commonly referred to as the KK masses m„: 



(54) 
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The interesting feature of (53) is that it separates into a four-dimensional part and a z dependent part. A mode of 
the form xo = e*'^''^''xo,n will solve the equation of motion (53) provided that the dispersion relation 



k^Ji^ + = 0, 



is satisfied, and hence modes labeled by n behave as four-dimensional massive particles. Technically, the advantage 
of working with Pq is that its (Euclidean) eigenvalues \n.k — k^k'^ + separate as a sum of a four-dimensional part 
plus the eigenvalue of the Schrodinger problem in the fifth direction. 

In the following subsections, we shall discuss how detPo is related to the determinant of our interest, detP, or 
more generally to detPg. For completeness, and in order to illustrate practical methods for calculating the effective 
potential, we shall consider dimensional regularization and zeta function regularization. Both methods will lead to 
identical results. 



C. Dimensional regularization 



A naive reduction to flat four- dimensional space suggests that the effective potential can be obtained as a sum over 
the KK tower: 



^2 



(55) 



Here D = 4 + l — eis the dimension of spacetime, and we have added (— e) dimensions parallel to the brane. The 
renormalized effective potential should then be given by an expression of the form 



rdiv 



(56) 



and the question is what to use for the divergent subtraction V^™ . Since Eq. (55) is similar to an ordinary effective 
potential in 4-dimensional flat space^, one might imagine that V can be obtained from just by dropping the pole 
term, proportional to 1/e; but this is not true for warped compactifications 

V{^i) i^V^ - (pole term). 

The point is that the theory is flve dimensional and the spacetime is curved, and this fact must be taken into account 

in the process of rcnormalization. 

Rather than proceeding licuristically from (56), we must take the definition of the effective potential Eq. (48) as 
our starting point, where it is understood that the formally divergent trace must be regularized and renormalized. 
In order to identify the divergent quantity to be subtracted, we shall use standard heat kernel expansion techniques. 
Let us introduce the dimensionally regularized expressions [14] 



J 



where 



Ce(s,£>)=1V 



2/z 



2s 



T{s) 



f 

Jo 



(57) 



(58) 



It should be noted that the operator Pg is positive and therefore the integral is well behaved at large ^. 

As is well known, the regularized potential Vf contains a pole divergence in the limit — > 5. To see that this is 
the case, one introduces the asymptotic expansion of the trace for small ^ [15,19,20], 



Tr 



OO 



(59) 



n=0 



^It should be mentioned also that each KK contribution in Eq.(55) is not just Uke a flat space contribution, because in warped 
compactifications the KK masses m„(ip, D) depend on the number of external dimensions parallel to the brane. 
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where a^^2 ^^"^ ^^'^ so-called generalized Seeley-De Witt coefficients. In (59) wc have introduced the arbitrary smearing 
function f{x). This is unnecessary for the present discussion, but it will be useful later on. For n < 5 their explicit 
form is known for a wide class of covariant operators, which includes our Pg. They are finite and can be constructed 
from local invariants (terms constructed from the metric, the mass term Ee and the smearing function /), integrated 
over spacetime. For even n, they receive contributions from the bulk and from the branes, whereas for odd n they 
are made out of invariants on the boundary branes only. 

For definiteness, let us focus on the simplest case of a Dirichlet scalar field, satisfying ^ 

X{z±) = 0. (60) 

We can use the result found in [16-18] to compute the Seeley-De Witt coefiicients for a Dirichlet field with a bulk 
operator P = — (□ + £). The lowest order ones for odd n are given by: 

«f/2(/,^) = ^^^4^E / vg'imd^-'x, (61) 



«?/2(/>^) = ^^384 ' ^id''-'x{f(96E-16n + 8nyy + 7)C^ - lO/C^./C'*'') + 0(/;„, (62) 

The most relevant for our purposes will be a^2- 



+ 15(7/C2 - lO/C^^/C'' + d.-Ryy - l&R) E 



+ ^Q-R? - 48 DTe - nnly - STlab'R-'''' + SlZabcd'^-'''""^ 

21 5 



47 215 



777 327 \ 1 

+ ^(^^"^^M.)' - —K.^.K.'-PK.p^K^^^) + Oif.y, /;,,,,)}. (63) 

Our notation is as follows, i;^ is a general scalar hmction. 'R-\^.,i = +^%c d ^ ' ' ' ^^'^ Ricmann tensor, TZbc = '^■"'bac 
the Ricci tensor and TZ = TZabg"''' is the curvature scalar. The extrinsic curvature is given by /C^j/ = {\/2)dyg^,y, where 
Qixviy) is the induced metric on y-constant hypersurfaces, and K, = IC^vg^"- The vector normal to the boundary 
is dy so the normal components are simply the y components. The a, 6, • • • indices run over the extra coordinate, 
and over the directions tangential to the branes, ^l, v, ■ ■ ■. The omitted terms, represented by 0{f.y, ...), are linear 
combinations of the derivatives of / with coefficients which depend on /C^^, E and its derivatives. 

As mentioned above, the integral (58) is well behaved for large ^. For small ^, the integral is convergent for 2s > D, 
as can be seen from the asymptotic expansion (59). In the end, we have to consider the limit s 0, and so we must 
keep track of divergences which may arise in this limit. For this purpose, it is convenient to separate the integral into 
a small ^ region, with ^ < A, and a large ^ region with ^ > A, where A is some arbitrary cut-off. Substituting (59) 
into (58), we can explicitly perform the integration in the small ^ region for 2s > D. This gives 



In §111 wc considered an even parity scalar field, which obeys Neumann boundary conditions. In fact, in the conformally 
trivial case we would get the same result for an odd parity field since the spectrum in the flat space problem is the same. The 
only difference is the presence of a zero mode in the even parity case, which does not contribute to the result (34) . 
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as,D) 



r{s) 



E: 

^n=0 



^n-D+2s 



(64) 



where we have used the standard notation 



The second term in curly brackets is perfectly finite for all values of s. Analytically continuing and taking the 
derivative with respect to s at s = 0, we have 



C'(0, D)^Y1 -^^^''nMPe) + finite, 

n=0 



(65) 



where the last term is just twice the integral in (64) evaluated at s = 0. Introducing the regulator e = 5 — D, the 
ultraviolet divergent part of Vg' is thus given by 



(66) 



The divergence is removed by renormalizing the couplings in front of the invariants which make up the coefficient 
a^2> so this infinite term can be dropped. The renormalized effective potential of our interest is therefore given 
by 



Ve = Um [ye° - 



divl 



(67) 



To proceed, we need to calculate Vg^, which in principle requires calculating a trace which involves the eigenvalues of 
Pg, and as mentioned above, these are not related in any simple way to the KK masses. 

However, it turns out that the dimensionally regularized Vg^ is independent of 9 when D is not an integer. The 
dependence of on 6 can be found in the following way. First we note that 



d0Ti 

where we have introduced 



Tr 



fe = delnQe, 



-^d^Tr fe{x)e 



(68) 



and the cyclic property of the trace was used. The above relation enables us to find the dependence of Vg' on the 
conformal factor: 



de \im dsCe{s,D) = \im ds—- / 

s^o s^o T(s) Jq 



(69) 



As with the expansion (64) we may again introduce the regulator A and separate the integral into a large ^ part with 
^ > A, which is finite and a small ^ part with ^ < A which contains the divergent ultraviolet behaviour. Assuming 
that 2s > D and integrating by parts, the resulting integrals in the small ^ region can be performed explicitly and we 

have 



de limdsCe{s,D) ~ lini ds . 

s— >0 s— >0 i [s) 



E;73:^T^<2(^,J'e) + finite 



n=0 



(70) 



As before, the last term just indicates the integral in the large ^ region. Provided that D is not an integer, all terms 
in square brackets remain finite at small s, and so the right hand side of (70) vanishes. Hence, we find that 



deVi" = 0, [D^ integer). 



(71) 



In other words, the dimensionally regularized determinant of Pg coincides with the dimensionally regularized deter- 
minant of Pq, and we have 



'P + ml{ipi,D) 



{D ^ integer). 



(72) 
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As was anticipated, we find that In Je vanishes in the dimensional regularization presented here in the sense given in 
(50). 

Finally, from (67) and (72), the renormalized effective potential is given by 

(73) 

where the Seeley-De Witt coefficient given in (63) with / = 1. The above equation bears the ambiguity in the 

choice of integration measure in the second term in square brackets. Different values of 6 give different results. If we 
take Qab as the preferred metric, then we should use 9 = 1, whereas if we take g^^^ as the preferred metric, we should 
use 6 = As we shall see in the next subsection, when we set D = 5 the coefScient 05/2 (Pg) is also independent 

of 9. Hence, the pole term in the second term in (73) is independent of 9, as it should, in order to cancel the pole in 
V^. However, the finite part does depend on the choice of 9. 

The right hand side of (73) is ready for explicit evaluation, which is deferred to the next section. We shall now turn 
our attention to the equivalent method of zeta function regularization. 

D. Zeta function regularization 

The method of zeta function regularization exploits the fact that the formal expression for the effective potential 

(57) is finite if the limit D ^ 5 is taken before the limit s 0. This can be seen from Eq. (64), where the term with 
n = 5 is finite if we set D = 5 before taking the derivative with respect to s and setting s 0. Clearly, the change in 
the order of the limits simply removes the divergent term V'^'^'" given in (66) and it reproduces the results obtained by 
the method of dimensional regularization (up to finite renormalization terms which are proportional to the geometric 
invariant o^'^^(Pe))- 

In zeta function regularization we define 

^0 = -^ ln(det Pe)^-^ lim d^e (s) , (74) 

where Cs{s) = C,g{s,D = 5) [see Eq. (58)]. As in the case of dimensional regularization, it is more convenient to 
calculate Vb than since the eigenvalues of Pq are related to the spectrum of KK masses. An important difference 
with dimensional regularization is that 

-2AdgVe = deCeiO) = 2a^,2{fe.Pe) ^ 0, 

a result which we already encountered in Ref. [7] (see also [16]). This can be seen from (70). If we set D = 5 from 
the very beginning, the term with n = 5 in Eq. (70) is linear in s, and its derivative with respect to s does not vanish 
in the limit s ^ 0. Here, and in what follows, we use the notation 

an/2 = «n/2- 

Integrating along the conformal path parameterized by 9, we can relate the effective potential per unit comoving 
volume Vo, with the "flat space" effective potential Vq as 

^^ = ^o-^f "5/2(/e'. ^e')- (75) 

The general expression for a5/2{fe, Pe) which applies to our case has been derived by Kirsten [16]. In Ref. [7] we 

evaluated the integral in (75) for the Randall-Simdrum case, in order to obtain Vi from Vq. Here we shall present an 
alternative expression for this integral which does not require the knowledge of a^/2{fej Pe)i but only the knowledge 
of a^i^{Pg) for dimension £) = 5 — e. This will also illustrate the relation between the method of zeta function 
regularization and the method of dimensional regularization. 

From the asymptotic expansion of the first and the last terms in Eq. (68), we have [21] 

dea^/2{Pe) = [D - n)a^,^{fe, P9). (76) 

Integrating over 9, we get 



Ve((f) = lim 

D^5 



(D-5) A "5/^ 



a?/o(-Pe) 
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{D - 5) /' a^j^{fe',Pe')de' = af,^{Pe) - af/^{Po). (77) 
Jo 

Writing D = 5 — e, we have 

'^5/2(/e'> ^e')'^^' = ^ [a^MPe) - <2(^o)] . (78) 

Note that, from (66) and (73), the previous equation can also be written as 

= Ve + V0^''' = Vo + Vf''. (79) 

This equation simply expresses the fact that the dimensionally regularized is independent of the conformal 

parameter 9, as we had shown in the previous Subsection [see e.g. Eq. (71)]. 

From (76), with D = n = 5, one finds that the coefficient a^/2{Pd) is conformally invariant [21], and therefore 

a^MPe) = a^/2{Po)- (80) 

Substituting this into (78), we obtain 

d 



f 

Jo 



a,/2ifo',Pe')de' = J^af/,{Pe) - ^^2(^0) (81) 



Thus, the integral in (75) can be evaluated in two different ways. One is by using the explicit expression of 05/2 (/, P) 
given by Kirstcn [16]. The other is by taking the derivative of the coefficients a^^^{Pg), given in (63) with / = 1, 
with respect to the dimension. [Note that the terms which are linear in derivatives of /, which we have just indicated 
symbolically in (63), disappear when / is a constant]. 

The quantity discussed above is nothing but the 9 dependence induced by the choice of integration measure. In the 
sense of (50), we have 

1 

lnJe = —J^ d9'a5/2{^,Pe'), (82) 

where we have used fg = de Inflg = (1 — 6')(/)/3c. Clearly, the effect of this factor is just adding local terms expressed 
solely in terms of (p and the metric to the classical action. The dependence of these terms is different from the change 
which results from a rescaling of the renormalization parameter fi. This corresponds to a shift in the coefficient of 
local terms proportional to 05/2 (-P)- 

V. EXPLICIT EVALUATION 

For simplicity we shall restrict attention to the case of massless fields with arbitrary coupling to the curvature: 

E = -^Tlg, 

and with Dirichlet boundary conditions. Here we shall use the method of dimensional regularization. Zeta function 
regularization is discussed in the Appendix. 
The eigenmodes in (54) arc given by 

Xn = Z^^'^{AiJ^{mnZ) +A2Y^{mnZ)). 

The index of the Bessel functions is given by 

u{D) = 1^1 - - m{D - 2)13 - 2]{i - UD)), (83) 



where 
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is the conformal coupling in dimension D. Imposing the boundary conditions (60) on both branes, we obtain the 
equation that defines implicitly the discrete spectrum of m„, 



F(m„) = J,,{m„r])Y^{rrin) - Y^{m„r])J^{m„) = 0, 



where we have defined 



V = 



(84) 



(85) 



The zeros of F are all real if p is real and rj is positive. Since r] is positive, the reality condition of v guarantees that 
all the KK masses are real. This provides a constraint for the possible values of ^ depending on q, 



16q{2 - 5g)' 
-(l-4g)^ 
16g(2 - 5q) ' 



q < 2/5, 
q > 2/5, 



(86) 



where we have used D = 5. Note that the values of ^ comprised between the minimal and the conformal coupling are 
allowed for any value of q. 




-0.25 




q=2/5 



FIG. 1. The shaded region corresponds to the allowed values in the {q -^) plane for a Dirichlet massless scalar field according to Eq. 
(86). Note that the range of ^ comprised between the minimal and the conformal coupling is allowed for any value of q 



In section V we concluded that the renormalized expression for the effective potential is 

1 1 



Veiifi) = lim 

D^5 



Consider first , given in Eq. (55). Performing the momentum integrations, we obtain 



where we have defined [7] 



y = + e/me - 4), 



C» =Y,m-^ = ^.jr'-^\nF{t) dt. 



(87) 



(88) 



(89) 



This regularized expression for the effective potential is finite when the real part of e is suflSciently large. In the last 
equation we have used that F{t) has only simple zeros which are along the real axis. The closed contour of integration 
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C runs along the imaginary axis, from t = +ioo to t = —ioo, skipping the origin through an infinitesimal path which 
crosses the positive real axis, and the contour is closed at infinity also through positive real infinity. 

Now the problem reduces to the computation of C, which can be done in the same way as in the case discussed in 
[7]. Skipping the detailed derivation, we simply give the final result: 

C(-4 + e) = -2/34(1 + 7?-^+^) - 2e(r?T)-2 {Ik + I/t* + t^V{t)) + 0{e^), (90) 

where 

/?4 = Y^(13-56i/2 + l6i/^). (91) 

Here we have introduced 

^^£<^ri/77, for g<l, 

z> Xn, for g > 1, 

to express the result for both g > 1 and q <1 cases simultaneously, where and z<_ are the largest and the smallest 
of Zj^ and Z-, respectively. Note that r < 1. 




FIG. 2. For q < I, in conformal coordinate, the singularity sits at z = 0, and the coordinate of the negative tension brane is smaller 
than the coordinate for the positive tension one, rj = zj^/z- > 1. For g > 1, the singularity is placed at z ^ oo and Z- > z+, so ?7 < 1. 



The constant coefficients Xk{i^), IiW) ar<^ calculable in principle, although their precise value is perhaps not very 
interesting since as we shall see these coefficients can be reabsorbed by finite renormalization. Their definition is given 
in the Appendix. Finally, V(t) is defined by 



V(t) = / 
Jo 



dpp^ln 1 



Kuirp) h{p) 



This provides the following expression for 

1 



V 



D 



(47r)2 



I 



K 



It 



V(r) 



+ 0(e), 



(92) 

(93) 
(94) 



where we have introduced the conventions Pi = (3i{D = 5) and (3'^ = d(3i{D) / dD\D=n-, since (3i depends on the 
dimension through (83). 

The next step is to subtract the divergent contribution. The Seeley-De Witt coefficient a^^^{Pg) can be computed 
from (63) for a generic dimension D, so we can expand the second term in the r.h.s of Eq. (73) as 
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{D - 5) A"^'^ 



(47r)2 



-Pi 



(95) 



where 5{6) is a constant whose precise value will be unimportant. 

The divergent parts of the two terms in Eq. (73) cancel, and the finite result for the effective potential per unit 
comoving volume is given by * 



Ve = 



(47r)2 



K 



V(r) 



Hpe + 1) 

(47r)2 



In 



■In 



(96) 



Pi 
'(47r)2 



ln(/i2;o)- 



Here, wc have elliminated some terms through redefinition of /i. 

As in the case of conformal fields, we must also allow for finite renormalization of the couplings in front of the 



invariants which make up the coefficient 05/2- These have the general dependence of the form aj_/C: 
these additions, we finally obtain: 



Ve{z+,z-) = 



Piipe + 1) 

(47r)2 



In {1J.1Z+) ln{iJ.2Z-) 

A I A 



+ 



(4. 



- [ 
Jo 



dx In 



/y(a;z<) Ki,{xzy) 



Ivixz^) K„{xz<) 



With 



(97) 



where fix and 112 are renormalization constants, and /3 and /34 are given in equations (31) and (91). This is the main 
result of this Section. In the limit of small separation between the brancs, (1 — r) <g; 1, the integral V behaves like 
(1 — t)~** (see Appendix C), and the logarithmic terms can be neglected. In this limit, the potential behaves like the 
one for the conformally coupled case, given in (33): 



Ve{z+,z-) ■: 



\z+ — Z- 



(98) 



For T <C 1, when the branes are well separated, the integral V behaves like t"^" and becomes negligible in the limit of 
small T (except in the special case when v is very close to 0). In this case, we have 



V0{z+,z_) 



Pi{pe + i) 

(47r)2 



ln(;m.g+) IjoAjj^Z-) 

A I A 





' / \ 21'" 


+ 


.(1) . 







(99) 



Due to the presence of the logarithmic terms, it is in principle possible to adjust the parameters fxi, ^2 so that there 
are convenient extrema for the moduli z+ and Z-. We shall comment on this point in the concluding section. 



VI. AN EXAMPLE 



Here we shall consider an example of a contribution to the moduli effective potential in the model of Lukas et al. 
[1,22,23], which may be relevant for the Ekpyrotic universe [2]. As mentioned above, this corresponds to the case 
q — 1/6. In principle, the fields in the action for this model do not have the canonical form, since in addition to the 
coupling to the metric they have unusual couplings to the dilaton (f>. Nevertheless, they can be studied along the 
lines of the previous Sections. For instance, the Heterotic M-theory model of [1,22,23] contains a scalar ^{x") in the 
universal hypermultiplet whose vev is zero in the background solution, and whose action is given by 

= - / y/gd'^x ^e-^{d^f. (100) 



''For 6» = 1, we find S = (-/3/3072){432-1936/3-900/32 +2257/3^ -1152(2+7(-l+/3)/3(l+3/3))$+3072/3(-2+3/3)(-2+7/3)^^}. 

*We can check that this coincides with Eq. (34) in the conformal case, f = 1/2 both for the Dirichlet or Neumann boundary 
conditions. Tlic trace anomaly of conformally coupled fields in this conformally flat background is zero for this topology and 
for any kind of boundary conditions, and so 04 = 0. Moreover, if f = 1/2 then Ki^2 = VW^pe"'' and I1/2 = yz/rrp sinhp. 
This means that Ik, 1i and V(r}) can be computed analytically, easily recovering the result Eq. (34). 
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This contains a kinetic term only, but it has a nonminimal form. However, changing to a new variable. 



we rewrite the action (100) as 

= ~ J \/—gd^x ^^dif^ — i (□(/) — ^ {d(j))^) if^^ + boundary terms 



(101) 



where the indices arc contracted with the five dimensional gab metric, and wc have integrated by parts in the first 
equality. For the Dirichlet case, the boundary terms 'generated' are not relevant since still the field (f vanishes there. 
The last equation shows that the potential term present in terms of tp mimics a nonminimal coupling to the curvature 
—^TZif^. Since the Ricci scalar for this background is 

^ 7 1 
9 t/2 

we conclude that the equivalent effective nonminimal coupling of is ^ = 15/28. Note that this point lies in the 
allowed region of values in the {q -^) plane defined by Eq. (86), which corresponds to an index of the Bessel functions 
1/ = 4/5. 

Hence, the contribution to the moduli effective potential induced by the field \1/ is given by (97) with v = 4/5. In 
this case we have f3 = 1/5 and /34 = —10179/80000. There are, of course, many other contributions, corresponding 
to all bosonic and fermionic degrees of freedom, ft turns out, however, that all bosonic contributions take a form 
similar to the one of the field \1/ [24]. If supersymmetry is unbroken, then these contributions are cancelled by the 
contributions from the fermionic degrees of freedom. But if the degeneracy between bosons and fermions is broken a 
la Scherk-Schwartz, for instance, then we expect that the resulting effective potential will be qualitatively similar to 
the one given by (97). A detailed study of this model is left for future research. 



VII. MODULI STABILIZATION 



In the limit of large interbrane separation, the potential (97) assumes a "Coleman- Weinberg" form for each one of 
the moduli, 



V{y+,y-) 



aK\yi)\n 



K{yi) 



(102) 



Here, we have expressed the potential in terms of the " curvature scale" 



so that K^{yi) behaves like a generic geometric operator of dimension 4 on the brane [such as the fourth power of 
the extrinsic curvature, or any of the operators in the integrand of Eq. (63)]. Working with instead of Zi has the 
advantage of relating directly to physical quantities, and hence it is easier to control whether we are in the range 
where the effective theory should be trusted or not. In particular, we should not allow Ki to be bigger than the cut-off 
scale of the theory. The constant a in (102) is given by 



(1 - 0)'l - 1 
(47r)2 



(103) 



where we sum over the contributions from all bulk fields %• The numerical coefficients (3^^ are given by Eq. (91). 

The value of 9 depends on the choice of integration measure in the path integral which defines the effective potential 
(see Section IV). If we adopt the measure associated to the Einstein frame metric gab which enters our original action 
functional (1), then we should take 6 = 1. However, this is not the only possible choice, as we have repeatedly 
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emphasized. The clasical action has a seahng symmetry whieh transforms both gab and the background scalar field 
(j). Using a conformal transformation which involves the scalar field, we may construct a new metric g)^fl which does 
not transform under scaling. If we adopt the measure which corresponds to this new metric, then we should take 
6 = \ — l/q. With this particular choice of 6 the coefficient a vanishes and the logarithmic terms in (102) disappear. 
Nevertheless, it is far from clear that this is indeed a preferred choice Here we take the attitude that the parameter 
6 is unknown, and that it should be fixed by a more fundamental theory of which (1) is just a low energy limit. 

The renormalization constants fii in (102) can be estimated by looking at the " renormalized coefficient" of the 
geometric terms of dimension 4 on the brane Ci{K) = ahi{K/ Hi). In the absence of fine-tuning, the Ci{K) are 
expected to be of order one near the cut-off scale K M, where is basically the five-dimensional Newton's 
constant. Hence, we expect 

~ Me-^*/«, (104) 

where Ci = Cj(M) ~ 1. In (102), we have also allowed for finite renormalization of local operators on each one of the 
branes. These operators are collectively denoted by 6ai. In order to ensure that the effective potential V does not 
severely distort the background solution, this correction to the brane tension must be much smaller than the effective 
tension of the brane in the classical background solution. Prom the Darmois-Israel matching conditions, this effective 
tension is of order M^Ki. Hence we require 

6ai < M^Ki < (105) 

In Section IV we considered contributions to the effective potential from massless bulk fields. These may have an 
arbitrary coupling to the curvature scalar of the standard form SJi-x^, or certain couplings to the background scalar 
field, such as the ones occurring in the Heterotic M-theory model considered in Section VI. However, if the model 
contains massive bulk fields, of mass m, then we expect terms proportional to rn^K^ in the effective potential. Even 
without massive bulk fields, we may expect the presence of lower dimensional worldsheet operators of the form M^K, 
M-^K^ and MK^, due to cubic, quadratic and linear divergences in the effective theory. Hence, we may expect that 
5ai has an expansion of the form 

6(Ti{Ki) ~ + juM^Ki + j2iM^Kf + ^siMKf + j^Kf + 0{Kf), (106) 

where Ki <C M, At <^ M and ju <C 1 in order to satisfy (105). For completeness, the above expansion includes the 
term proportional to Kf. It should be understood that this term is only present in the particular case 9 = 1 — 1/q 
(corresponding to the scale invariant metric), since for other values of we assume that it is reabsorbed in a redefinition 

oi ^i^ [see Eq. (104)]. 

The local terms may in principle stabilize the moduli at convenient locations. Note that this effect is due to the 
warp factor and vanishes in flat space (where the coeflicients /34 vanish). The effect also vanishes accidentally in the 
RS case, because the curvature scale K{y) is constant. The position of the minima are determined by dy^V = 0. This 
leads to the conditions 



5a, = -Kf 



(i-,)i„if).i 



+ K,'-^. (107) 



where the prime on Sa, indicates derivative with respect to K,. Also, we must require that the minima occur at 
an acceptable value of the effective cosmological constant. Using the condition (107), we can write the value of the 
potential at the minimum as 

Vmin = §E (^) +a + Kr^6a'^ < IQ-^^X- (108) 

The latter condition will require one flne tuning amongst the parameters in (106). 



^In particulax, as mentioned in Section II, for certain values of the model parameters our action can be obtained by dimensional 
reduction of a 5-l-n dimensional Einstein-Hilbert action with a cosmological term. Since in the higher dimensional theory only 
gravity is present, there is only one possible choice for the metric and the issue of the measure docs not arise. Preliminary 
analysis of the case with n=l [24] indicates that the logarithmic terms do indeed arise in the limit when the radius of the sixth 
dimension is much smaller than the interbrane separation. For this case the relevant value of the parameter is 6 = 0. 
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An interesting question is whether the effective potential (102) can generate a large hierarchy and at the same time 
give sizable masses to the moduli. As discussed in Section II. B, the hierarchy is given by 



m 

mrii 



M 



K- 



2q-2 



(109) 



where m ~ TeV is the mass of the particles which live on the negative tension brane, as perceived by the observers on 
the positive tension brane. Consistency with Newton's law at short distances requires K+ > {TeV)^ /rrip ^ (mm)^^, 
and consistency of perturbative analysis requires K_ < M. With these constraints, the observed hierarchy h ^ 
exp(— 37) can only be accomodated for > 5/4. To proceed, we shall distinguish two different cases. 



A. Case a: 



This is the generic case, where the coefficients 71^, 72^ and 7.3^ in the expansion of 5ai{K) are not too suppressed. 
In this case, the logarithmic terms in the effective potential are in fact subdominant, and the minima of the effective 
potential are determined by AqSai ~ Ki8a\. 

The present discussion applies also to the special case where 9 = 1 — 1/q (corresponding to the measure associated 

with the scale invariant metric S'^t''), so that no logarithmic terms are present in the effective potential. Note that 
terms of the form 'juKf and ^uM^Ki in the expansion of 5ai{K) [see (106)] cannot be avoided. The first is necessary 
in order to renormalize the effective potential, and the second is already present at the tree level, so it just corresponds 
to a shift in the existing parameters in the classical action. 

Quite generically, this will lead to stabilization of the moduli near (or slightly below) the cut-off scale Ki = XiM, 
with Xi 1. Hence we have 

/i^ ~ exp[2(g- l)ln(A+/A_)]. 

Since the logarithm is of order one, an acceptable hierarchy can be generated provided that q > 10. This is "close" to 
the RS limit q 00. In this case. m.p ~ M. On the positive tension brane the parameter A+ has to be fine tuned so 
that the effective cosmological constant is 122 orders of magnitude smaller than the Planck scale. A straightforward 
calculation shows that the physical mass eigenvalues for the moduli ip± in the present case are given by 

2 —2 —2t^4 <^ 2 2 — li,2 — 2 r^4 <' i,2 2 

m+ ~ q nip ^ nip , ml ~ q h nip K_ ^ h nip. 
Thus, the lightest radion has a mass comparable to the TeV scale. 



B. Case b: 



This corresponds to the case where almost all of the operators in (106) are either extremely suppressed or completely 
absent, due perhaps to some symmetry. In particular, we shall concentrate on the possibility that 

since an operator of this form is already present in the classical action (1), and it is the only one in the expansion 
(106) which is allowed by the scaling symmetry. In this case, and assuming for simplicity that the negative tension 
brane is near the cut-off scale K- M, we can rewrite (108) as 

^min - { (in (i^+^+) + /.«(-^)/(^-^) (in (K_^_) + i) I 

Here we are assuming that 6 ^ 1 — 1/ q {so that a ^0), since the alternative case was already discussed in the previous 
subsection. For q > 1, the first term dominates and the condition of a nearly vanishing cosmological constant forces 
« /i_)_e~^/^. A fine-tuning of A+ will be necessary in order to satisfy the condition (107) for such value of K^. 
The hierarchy is given by 
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where /i+ is given by (104). Since the effective coupUng a can be rather small, a large hierarchy may be obtained 
even for moderate g > 1. A straightforward calculation shows that at the minima of the effective potential (102) 
d^^V = 12a(l + 2q)-^aXKlip+\ and dl_V ^ aq ^a^K^ipJ^ . Hence, we find that the physical masses for the 
moduli fields and which appear in (15) are given by 

ml ~ ag-2/ii2/(2«-i) ^2 ^ ^ c.g-i/i2+V(2«-i)m^. 

Associated with the eigenvalue m+ there is a Brans-Dicke (BD) field with BD parameter lubd = — 3(j'/(l + 2q). 
Therefore, we must require m+ > {mm)~^, which in turn requires g > 2. A stronger constraint on q comes from the 
eigenvalue m_ , since the corresponding field is coupled to ordinary matter with TeV strength. The mass of this field 
cannot be too far below the TeV, otherwise it would have been seen in accelerators. This requires q to be rather large 
9^ 10. 

VIII. CONCLUSIONS 

We have studied a class of warped brane- world compactifications, with a power law warp factor of the form 
a{y) = {y/yaY '^^'^ ^ dilaton with profile cj) oc ln(?y/?yo)- Here y is the proper distance in the extra dimension. In 
general, there are two different moduli y± corresponding to the location of the branes. (in the RS limit, g — » oo, a 
combination of these moduli becomes pure gauge). 

Classically, the moduli are massless, but they develop an effective potential at one loop. We have presented methods 
for calculating this effective potential, using both zeta function and in dimensional regularization. The procedure used 
in the case of zeta function regularization is a simpler version of the one we introduced in Ref . [7] . Our treatment of 
dimensional regularization formalizes and extends the approach adopted in Refs. [8-10]. An important point is that 
the divergent term to be subtracted from the dimensionally regularized effective potential is proportional to the the 
Seeley-de Witt coefficient 05/2, given in (63). In the RS model, this coefficient behaves much like a renormalization 
of the brane tension, but it behaves very differently in the general case. 

In general, the effective potential induced by massless bulk fields with arbitrary curvature coupling is given by (97). 
In the limit when the branes are very close to each other, it behaves like V cx a'*|?y-|- — y-l"**, corresponding to the usual 
Casimir interaction in flat space. Perhaps more interesting is the moduli dependence due to local operators induced 
on the branes, which are the dominant terms in V{y+,y-) when the branes are widely separated. Such operators 
break a scaling symmetry of the classical action, which we discussed in Section II, but nevertheless are needed in 
order to cancel the divergences in the effective potential. If we denote by K{yi) — q/yi the extrinsic curvature of the 
brane at the location y — yi {i = ±), a renormalization of the brane tension parameters a± in the classical action (1) 
induces terms proportional to a{yiYKi in the effective potential. These terms scale like the rest of the classical action 
under the global transformation (19-20). On the other hand, the divergences in the effective potential, proportional 
to the coefficient a;,/2{P), require world-sheet counterterms which arc proportional to a{yi)^K^{yi). These have the 
wrong scaling behaviour [they simply do not change under (19-20)] and hence they act as stabilizers for the moduli. 

In addition, there are terms proportional to a{yiY K^{yi)<p{yi). The coefficient in front of the latter terms depends 
on the choice of the measure in the path integral. Different choices are possible, which are related amongst each other 
by dilaton dependent conformal transformations. Because of the conformal anomaly, different choices are inequivalent, 
but they are simply related by the addition of local world-sheet operators to the action. These are given by the r.h.s. 
of (82). Since the scaling symmetry of the action is anomalous, the presence of these local terms could perhaps have 
been guessed without the need of detailed calculation. Since ((> behaves logarithmically, these terms have the form of 
Coleman- Weinberg type potentials for the moduli yi, and they can also act as stabilizers for the moduli. 

To conclude, we find that worldsheet operators induced on the brane at one loop easily stabilize the moduli in 
brane-world scenarios with warped compactifications, and give them sizable masses. If the warp factor is sufficiently 
steep, q > 10, then this stabilization naturally generates a large hierarchy, as in the Randall-Sundrum model. In 
this case, the mass of the lightest modulus is somewhat below the TeV scale. This feature is in common with the 
Goldberger and Wise mechanism [25] for the stabilization of the radion in the RS model. For q < 10, the stabilization 



Here we arc considering the situation where the mass of tp- is much larger than the mass of tp+, and where the visible 
matter is on the negative tension brane. In this case, since y- = const., visible matter is universally coupled to the metric g^i., 
and the BD parameter corresponding to ip+ can be read off from (15). Alternatively, it can be read off from (A4) after setting 
y- = const. 
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is also possible, but if wc also demand that the hierarchy /i ^ 10 is generated geometrically, then the resulting 
masses for the moduli would be too low. 
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APPENDIX A: ANOTHER FORM OF THE MODULI ACTION 

In this Appendix, wc present alternative forms of the moduli action in terms of the metric induced on either brane. 
This highlights the fact that the interactions with matter on the branes are of Brans-Dicke type. 

Changing variables in the background solution (2) according to y = y_ + RQ, with < O < tt, we can recast 
the moduli as the radius of the orbifold R and the 'would be' distance ?/_ from the negative tension brane to the 
singularity. Then the line element of our solution is 



yo 



rj^.^dx^dx", (Al) 



with R and y_ arbitrary constants. 

In order to take into account the gravitational degrees of freedom, wc promote R and y_ (and the induced 4d 
metric) to 4d 'moduli' fields depending on x^, so that the ansatz we shall consider takes the form 

ds^ = R{x^'f dO^ + (1 + rix^") sf'' g'^^'Ux^dx'' , (A2) 
<f> = V6q [In (1 + r{x^) 6) + In {v-{x^)ly^)\ (A3) 

Here we have absorbed y_(a;'*) into the 4d metric glj^ ■, and defined r = R/y-. We insert this ansatz back to the 
action (1) and integrate over © to get the kinetic terms of the 4d effective action. In terms of the metric on the 
negative tension brane we find 



2 



where indices are contracted with the ^ metric, TZ'^ ^ is the 4d curvature scalar related to it, and we defined 



Note that A and y_ have diagonal kinetic terms, with A depending on the hierarchy of scales between the branes only. 

The contribution to the kinetic part of the effective action from these moduli as seen from the positive tension brane 
is obtained from Eq. (A4) by writing gj^^ as A~^''.g^i^ (equivalently, we can change A by 1/A and y_ by y+ = Ay_). 
After partial integration and a global change of sign, wc get 



^6 



z^/..^^(i^)„., 

.3qX-(^^+^)-:-^l (A6) 




2q + l ) yl +^^^ A2 
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with indices contracted with g^^K 

It is easy to see that both A and y_ or j/+ can be put in the form of Brans-Dicke (BD) fields, with a generic action 

Ignoring the global A-dependent factor, y- and y+ are pure Brans-Dicke field up to a rescaling, with Brans-Dicke 
parameter u)^j-, = 3q. In the case of A, again ignoring the global y_ (or factor, the Brans-Dicke parameter are 

^Bi,^^(A) = ^(AW)-l). (A8) 

From the positive tension brane, and for g ~ 1, the observational limit, ujbd > 3000 [26], is fulfilled with A ~ 15. 
This requires an interbrane proper distance ttR at least about 1 order of magnitude larger than the distance from the 
singularity at y = to the negative tension brane, at y = In the negative tension brane, we find that the BD 
parameter is negative but greater than —3/2 for any A or q, as happened in [27]. We reproduce the results of [28] in 
the limit of g ^ oo. 



APPENDIX B: DIMENSIONAL REGULARIZATION VERSUS ZETA FUNCTION REGULARIZATION 



1. Zeta function regularization 



Let us rederive the main result of Section V using zeta function regularization. For simplicity, we shall restrict 
attention to the case 9 = 1. In this case, we have explicitly calculated the integral (78) along the conformal path 
using both methods described in Section IV. This exercise leads to the result 



V = Vo 



(47r)2 



1 1 



.4 ,4 



(Bl) 



where a = (/3/3072){144 + /3(784- 1692/J + 335/32 - 192(3/3 - 14) (3/? - 2)^)}. The zeta function associated with the 
operator Pq is given by 



d*k 
(2^ 



E 



Performing the momentum integrals in (B2), we have 

2.^2.-4p(^ _ 2) 



Co(s) = ^- 



(47r)2r(s) 



C(2s-4), 



Substituting (B3) into (74), we have 
Vo{z+,Z-) = 

and from (90) we obtain 



1 



2(47r)2^4 



ln(/.^_) + -K(-4) + C'(-4) 



Vq{z+,Z-) = 



1 



(4^) 



II , V(t) 



f34 



(47r)2 



Z I z 



(B2) 



(B3) 



(B4) 



(B5) 



Substituting in (Bl) we recover Eq. (97) up to finite renormalization of /z. 
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2. Vo in dimensional regularization 



We shall now reproduce the result for Vq by using dimensfonal regularization. Again, this is a redundant exercise: 
the calculation of an effective potential (be it V or Vo) will give the same answer whether it is done in dimensional 
or in zeta function regularization. Nevertheless, it is interesting to do it explicitly since this calculation is closest in 
spirit to the standard flat space calculations in four dimensions. 

Adding up the dimensionally regularized effective potential per comoving 4-volume due to all KK modes, we have 

Performing the momentum integration for each mode, we obtain 

where we used to„ defined in (85). This regularized expression for the effective potential is finite when the real part 

of e is sufficiently large. Performing analytical continuation in e, the summation over KK modes ^^"^ can be 
identified with the zeta function C(— 4 + e). The pole part proportional to 1/e is identified with 



Subtracting this divergent part, we get the renormalized expression for the effective potential as 

1 



Vo = V""^ - V^'^ = - lim 



2(47r)2 



(B9) 



Consequently, we find that the dimensional regularization method reproduces the previous result (B5). 

Here, one remark is in order. In a usual 4-dimensional problem, the divergent part is given by the Seeley-De Witt 
coefficient a^- In the present case, since the background is 4-dimensional flat space, consists of only one term 
proportional to m^^ for each KK mode. Subtraction of this counter term for each mode leads to the expression 
X]„TO^logm„, which is still divergent. This is not a surprising fact because the problem is essentially 5-dimensional 
as it is indicated by the existence of the infinite tower of the KK modes. Usually, this point is bypassed in the 
literature by evaluating the divergent sum with the help of the generalized zeta function. 



3. Alternative regularization 

In section V, we were interested in the class of conformally related operators Pe, and therefore it was important to 
do a dimensional extension of the spacetime such that all geometries labeled by 6 would be conformally flat. In other 
words, we added dimensions "parallel" to the brane, whose "size" was also affected by the warp factor. It should be 
stressed that this was done for computational convenience, since is independent of the parameter 6 only in this 
regularization. Putting aside computational considerations, nothing prevents us to extend the spacetime in any way 
we please, and the results should still be the same. To illustrate this point, let us consider an alternative dimensional 
extension of our 5-dimensional curved space M. to a, simple direct product space given by M^*^ x Ad. The dimensional 
regularization is done by an analytic continuation of the number of added dimensions e, while the manifold M. is kept 
unchanged. 

Since this is a direct product space, the eigenvalues of the (5 — e)-dimensional d'Alembertian are given by a simple 
summation of those in each space, M. and 

A(5_e) =kl + \. 

Then, the dimensionally regularized effective potential per unit comoving volume is given by 




(BIO) 
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This quantity is evaluated by introducing the function 



where Ci is the zeta function defined in (58) with 9 = 1. Then, using (Bll), we can rewrite (BIO) as 

= "2^^'^°^ = -2^(l/47r)-/2r(e/2)Ci(e/2). (B12) 

The above expression contains the object Ci, which wc have encountered in zeta function regularization. However, it 
should be noted that now the regularization parameter is not s, the argimicnt of the function T, but the dimension 
of the product space e. Therefore, even after we take the limit s ^ 0, V"^^ still diverges as 1/e. 

The divergent piece in the dimensional regularization in the D dimensional problem is the Seeley-De Witt coefficient 
a/3/2. Prom (58) and (59), it is easy to relate this divergent piece with the value of Ci(0) as^^ 

^^^'^ = ^^^ = ^Ci(0). (B13) 

Subtracting this divergent piece from (B12), we obtain the renormalized value 

V = - ydiv = -i^ (c; (0) - {In (l/47r) + 7}Ci(0)) , (B14) 

where 7 is Euler's gamma. This coincides with Eq. (74) for ^ = 1 up to a redefinition of fi (note that with our 
conventions, Cl(0) depends on /x.) 

APPENDIX C: DEFINITION OF THE CONSTANTS Ii AND Ik 

The result (90) involves the two constants 1/ and Tj. Although the actual value of this constants seems to be 
irrelevant, since they can be shifted by finite renormalization of worldsheet operators, we will nevertheless for the sake 
of completeness give their definition in terms of the parameters of the model. 

We can express the asymptotic expansion of the combinations of modified Bessel functions used in subsection IV as 

p °° Ai) 

Since both ^q^^ = 5q^^ = 1, and the obsevration that P^^^ = (3^^ = /J4, we can define coefficients fir for both K^, and 
Iv as 



Defining the functions 



,,^0 P J 7^1 P' 



3 f>(.K) o 



^^"np) = E^ + 3 



.=1 p' 



r=l 



'^^Strictly speaking aP^j has to be evaluated in the regularized (5+e)-dimensional space. However, since the added e dimensions 



^5/ 

are trivial, it is identical to 05/2(-Pi). 
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the coefficients Xj and Tk are defined as 



1i = -7/34 + ^ dpp" (in [^2^ e-^/,(p)] - , 
Tk = -7/34 + ^ rfpp' (in [^/2^e''K^{p)'\ - . 



(C4) 



APPENDIX D: ASYMPTOTIC FORM OF V(t) 

The behaviour of V(r) defined in (92) for r <C 1 is given by 



V(r) = - 



2v 



0{t^ In: 



(Dl) 



This corresponds to a large separation between branes. In this limit the integral is generically negligible compared 
with the logarithmic terms in (97). 

The limit of small separation between branes corresponds to 1 — r -C 1. In this limit, the integral can be approxi- 
mated by taking the arguments of the Bessel functions to be large. Using the asymptotic expansion 



K,{Tp)I,{p) 



(D2) 



we have 



j dpp^ In (l - 

- 1 r 



dx 



- 1 



3C(5) 
8(1 -t)4- 



(D3) 



In the second equality, we performed integration by parts and a change of variable. Here, C, is the usual Riemann's 
zeta function (^(5) = 1.03693 • • •. Using the relation 



(D4) 



and substituting in (97), we find that in the limit of small brane separation, the effective potential reduces to the one 
we had found in the massless conformally coupled case, given in equation (33): 



V{Z+,Z.)r.-- 

\z+ - z. 

Equation (D4) is a particular case of a more general formula 

(2n)! 



|4 ■ 



(D5) 



C'(2n) 



C(2n + 1), 



22n+1^2n 

valid for positive integer n. This can be derived from the perhaps better known relation [29] 

C(l -z) = 2i-^7r-^C(-2)r(^) cos(7r^/2), 
by setting 2; = 2n + 1 after differentiation of both sides of the equation with respect to z. 
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